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! Abstract 

CO 

\ A next-to-leading order analysis of inelastic electroproduction of charm is 

ON ' performed using an interpolating scheme which maps smoothly onto massless 

• QCD evolution at large Q 2 and photon-gluon fusion at small Q 2 . In contrast 

Oh, with earlier analyses, this scheme allows the inclusion of quark and target 

Oh. mass effects and heavy quark thresholds, as well as possible non-perturbative, 

r- 1 \ or intrinsic, charm contributions. We find no conclusive evidence in favor of 

an intrinsic charm component in the nucleon, although several data points 
which disagree with perturbative QCD expectations will need to be checked 
by future experiments. 



I. INTRODUCTION 

Understanding the role played in the nucleon by heavy quarks, such as the charm quark, 
is necessary for a number of reasons. Firstly, one cannot claim to have unraveled the rich 
structure of the nucleon sea until one has mapped out the details of the distribution of 
its virtual charm and heavier^] flavors. Secondly, in the absence of a direct probe of gluons, 
charm leptoproduction remains one of the main sources of information on the nucleon's gluon 
distribution. Furthermore, tagging charm in neutrino and antineutrino scattering allows one 
to probe the strange and antistrange quark densities in the nucleon. From a more theoretical 



1 Although in practice the direct accessibility of bottom and top quark densities is likely to remain 
elusive for some time. 
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point of view, in order to have a reliable procedure through which to analyze deep-inelastic 
scattering data, one needs to consistently incorporate heavy quark masses and threshold 
effects in the QCD evolution equations. 

Recently there have been important new data on the charm structure function, F%, of 
the proton from the HI and ZEUS Collaborations at HERA, which have probed the 
small-x region down to x = 8 x 10~ 4 and 2 x 10 -4 , respectively. At these values of x, the 
charm contribution to the total proton structure function, F% , is found to be around 25%, 
which is a considerably larger fraction than that found by the European Muon Collaboration 
at CERN H at somewhat larger x, where it was only ~ 1% of F%. Extensive theoretical 
analyses in recent years have generally served to confirm that the bulk of the F% data can 
be described through perturbative generation of charm within QCD. 

At the same time, there are several lingering pieces of data which seem to suggest the 
possibility that a small component of charm exists which is intrinsically non-perturbative 
in origin One of these is the EMC data at large x, some of which appear to 

lie above the perturbative QCD predictions, and which have in the past been taken ||] as 
evidence in favor of a non-perturbative, or "intrinsic" charm component. Furthermore, as 
recently discussed in Ref . |J , there are some indications of intrinsic charm also from hadronic 
reactions, such as leading charm production in nN and Y N scattering. As found in Ref. M , 
some intrinsic charm may account for the larger than expected number of fast correlated 
J/ijj pairs seen in the NA3 irN experiment at CERN ||10|| , as well as for the anomalous 
polarization of J/ ip seen in inclusive J/ ip production in ttN collisions [[11] . 

The initial analyses of the leptoproduction data suggested an intrinsic charm component 
of the nucleon with normalization of around 1%, although subsequent, more sophisticated, 
treatments incorporating higher-order effects tended to find only a fraction of a percent. 
The basic starting assumption in these analyses is that the perturbative generation of charm 
proceeds through the photon-gluon fusion process, illustrated in Figure 1. In their pioneering 
analysis, Hoffman and Moore [BJ considered the 0(a s ) corrections to the intrinsic charm 
distributions, as well as quark and target mass effects on the charm cross section. With 
these refinements, their work suggested that an intrinsic charm component of the proton, at 
a level of order 0.3%, was allowed — even required — by the data. In a more recent study, 
Harris, Smith and Vogt || reanalyzed the EMC data within the framework of Ref. || with 
updated parton distributions, and 0(a s ) corrections to the hard scattering cross section. 
They found that while no clear statements about intrinsic charm could be made at lower 
energy transfers, v < 100 GeV, data at larger v ~ 170 GeV could best be fitted with an 
additional 0.86 ± 0.60% intrinsic charm component ||. 

Despite these claims of evidence for intrinsic charm, a number of recent global parame- 
terizations of parton distributions [I3|.|R| have concluded that all _F| data, including those 
from the EMC, can be understood in terms of perturbative QCD alone, with no additional 
intrinsic charm component necessary. In view of these conflicting results, it seems impor- 
tant that the question of the presence or otherwise of intrinsic charm be addressed using 
the latest available theoretical techniques for calculating F£, as well as different models for 
intrinsic charm. 

To set the context for our analysis, we should point out a number of features which have 
been common to all earlier treatments of intrinsic charm in leptoproduction. Firstly, one 
invariably assumes that photon-gluon fusion (together with perturbative QCD corrections) 
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is solely responsible for perturbatively generated charm, irrespective of the scale Q 2 at which 
one works. In light of recent developments in incorporating quark mass effects in perturbative 
QCD, it is now in principle possible to go beyond this approximation and include more 
appropriate treatments in kinematic regions where photon-gluon fusion alone may not be 
the only relevant perturbative mechanism. A second point concerns the way that intrinsic 
charm is combined with the perturbative contributions. In particular, no treatment to date 
has considered the effects of Q 2 evolution of the non-perturbative charm distribution from 
the "intrinsic" scale to the scale where data are taken. 

In this paper we perform a comprehensive next-to-leading order analysis of the charm 
production data in which we compare several different methods of implementing charm in 
QCD structure functions. In particular, we consider: 

• Massless Evolution. This approach assumes that below a certain scale, /i 2 , there are rif 
active flavors. Above this scale, one switches the number of flavors to rif + 1 in the coupling 
constant, coefficient functions and in the parton distributions. This scheme is sometimes 
referred to as the Variable Flavor Number Scheme (VFNS). 

• Photon-Gluon Fusion. In this scheme, the heavy quark is never treated as a parton, 
in the sense that the number of active flavors does not change, and a quark distribution 
for the heavy quark is never introduced. The entire contribution to the physical structure 
function from heavy quarks is generated through the photon-gluon fusion process depicted 
in Figure 1. This scheme is usually referred to as the Fixed Flavor Number Scheme (FFNS). 

• Interpolating Schemes. While photon-gluon fusion should provide a reliable description 
of charm production in the region close to m 2 , it certainly should fail at large scales because 
of the lack of resummation of the logarithms in m 2 . Conversely, massless evolution should 
be good at large Q 2 , but not in the region around m 2 . Schemes such as that proposed by 
Aivazis et al. [0] reproduce the relevant features of both the VFNS and FFNS approaches 

- namely, a reduction to FFNS when Q ~ m 2 , and to VFNS when Q 2 > m 2 . Unlike Ref. 
fH]] , however, we include the full quark mass dependence not only in the gluon sector, but 
also in the quark sector, which has become possible since the recent work of Kretzer and 
Schienbein ||15|| . 

Throughout the analysis we shall work in the MS scheme, and introduce an explicit 
charm distribution at a scale m 2 for the VFNS and interpolating schemes. Furthermore, our 
philosophy is that if there is any intrinsic charm at all in the nucleon, it should be considered 
as a parton distribution for all scales Q 2 > m 2 , if one is to use MS anomalous dimensions 
when evolving the distributions. 

In the following sections we describe in turn the three schemes for incorporating charm, 
beginning with the standard massless QCD evolution and the photon-gluon fusion process 
in Sections II and III, respectively. The interpolating scheme connecting these two limiting 
cases is described in Section IV. Incorporation of intrinsic charm into the analysis is outlined 
in Section V, where we survey several models for the non-perturbative charm component. 
Finally, in Section VI some conclusions are drawn from the analysis. 
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II. MASSLESS EVOLUTION (VFNS) 



The simplest approach to describing charm in the nucleon is to assume that the charm 
distribution is c(x, /i 2 ) = c(x, /i 2 ) = below a certain scale, /i, and evolve charm quarks 
as massless partons. Indeed, until recently, this has been the standard approach adopted 
in global analyses of parton distribution functions. In a scheme such as MS, the scale /i 
is taken to be the charm quark mass, m c ||16||. Above threshold, Q 2 > m 2 , the number of 



active flavors increases by one in the coupling constant, splitting and coefficient functions. 
In this way the charm quark distribution is generated solely through radiative corrections 
from gluons and light quarksQ present at Q 2 < m 2 . 

At next-to-leading order, the charm quark contribution to the proton structure function 
is given by: 

F c 2 (VFNS) (x, Q 2 ) = e 2 C q {x, Q 2 ) ® (xc(x, Q 2 ) + xc(x, Q 2 )) + e 2 C g (x, Q 2 ) ® xg{x, Q 2 ) (1) 

with the boundary condition c(x,m 2 ) = c(x,m 2 ) = 0. The convolution "®" here is defined 
by f(x) (g) h(x) = /J dyf(x/y)h(y). The quark and gluon coefficient functions are those for 
massless quarks: 

C q (x, Q 2 ) = 6(1 -x) + ^f^C F L ( l ^—A\ _ - 2(1 + X ) ln(l - x) 

4tt { \ 1-x J + (l-.x)+ 

1 + x 2 1 

-2 Ins + 6 + 4a: L (2a) 

1 — x J 

C g (x, Q 2 ) = ^4^4T F ((1 - 2x + 2x 2 ) In - 4 + 32x(l - x)\ , (2b) 



where the (1 — x) + refers to the standard "-(-"prescription [|T^. 

The gluon coefficient function, C g , is a pure singlet quantity, and as such is proportional 
to the number of active flavors. Since one is calculating the contribution to from a 
single flavor, C g defined in Eq.(^b|) is the usual gluon coefficient function divided by the 
number of active flavors. The gluon distribution is evolved by taking as input a gluon 
distribution function at a scale m 2 , together with a quark singlet distribution, E(x, m 2 ) = 
J2 q =u,d,s,c (l( x y m l) + ?( x ? m c))- Using the boundary condition c(x, m 2 ) = c(x, m 2 ) = 0, the 
evolution is performed with 4 flavors. 

The evolution of the charm quark density is slightly more involved. Firstly, one evolves 
the singlet combination from m 2 to Q 2 using the full number of active flavors. The quark 
non-singlet combination, rj(x, Q 2 ) = Y, q = u ,d,s {q{ x -> Q 2 ) + <l( x , Q 2 )) ~ 3 (c(x, Q 2 ) + c(x, Q 2 )), 
is then evolved over the same range. Since c(x,m 2 ) = c(x,m 2 ) = 0, it turns out that one is 
in practice evolving J2 q=u ,d,s (<l(. x , m c) + Q.i x i m D) fi rs f as a singlet and then as a non-singlet. 
The charm distribution at Q 2 is then given by: 

c(x, Q 2 ) + c(x, Q 2 ) = - (Z(x, Q 2 ) - r)(x, Q 2 )) . (3) 



2 Just as for the c quark, the u, d, and s quarks are light, so too for heavier quarks, such as b, the 
u, d, s and c quarks would be light, and so on. 
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For large Q 2 , Q 2 3> m 2 , one expects the generation of charm through massless evolution 
to be a good approximation, a result confirmed in previous studies [18|. However, at smaller 

Q\ Q 2 



mass effects become important, and this approach eventually breaks down. 
In this region a better approximation is the photon-gluon fusion process, which we discuss 
next. 



III. PHOTON GLUON FUSION (FFNS) 

In the photon-gluon fusion approach the charm quark is treated not as a parton, but 
rather as a heavy quark, so that a charm quark distribution is never explicitly introduced. 
Hence, c = c = for all x and Q 2 . The charm structure function is given directly by a 
convolution of the gluon distribution and the hard photon-gluon cross section, Figure 1. 

At lowest order, the only process producing a heavy qq pair is the photon-gluon fusion, 
and the structure function is given by: 

,c (FFNS), n2 ^ _ 2 «s0 2 ) [ Zma * ^ u ^ 2 X _ ( X 



F^-»>(x,Q 2 )=e 2 c ^ ^ dz H g (z,Q 2 } m 2 c ) ~-g (j-, , (4) 

where z max = 1/(1 + 4m 2 /Q 2 ) and fi 2 is the mass factorization (and renormalization) scale. 
The partonic cross section for producing a massive quark pair is given by: 

HJz, Q 2 , m 2 c ) = 4 ( 1 - 2z + 2z 2 + 4z(l - 3s)^f - 8z 2 ^$) h/ 1 + ^ 



Q 2 Q*J \l-f3 / 
Q 2 



+ / 3{32z(l-z)-4-16z(l-z)-f , (5) 



with 



Az ml , , 



Note that the coupling constant is calculated at /i 2 with 3 light flavors, the same number 
as used throughout in H g . There is no evolution equation in this approach — the distri- 
bution g(x,fi 2 ) is never evolved. Evolution equations are introduced only when the logs in 
Q 2 /m 2 are resummed via the renormalization group equations. In that case a heavy quark 
distribution is introduced and one recovers Eq. ([[]). 

One should mention that the gluon distribution in Eq.([|) is that given by a leading- 
order analysis of the data. This is because in Eq.(^j) one has terms of the form 
a s (fj, 2 ) \n(Q 2 /m 2 ) g(x, /x 2 ), which are typical next-to-leading order contributions if the gluon 
distribution is calculated in leading order — hence one is calculating 0(a s ) corrections to 
the structure function. The use of a next-to-leading order gluon would produce next-to-next- 
to-leading order corrections, which would then require calculation of radiative corrections to 
the photon-gluon fusion term, and the addition of extra terms in Eq. (|4]) arising from heavy 
pair creation from light quarks [|T^]. 

The advantage of the photon-gluon fusion formulation is its simplicity. However, the 
use of Eq.(f|) could be problematic at very large Q 2 because of the large ln(Q 2 /m 2 ) term 
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appearing in H g . One way of circumventing this difficulty is to construct a scheme which 
retains the features of photon-gluon fusion at moderate Q 2 , but maps onto massless evolution 
when Q 2 becomes very large. 



IV. INTERPOLATING SCHEME 



The description of charm production from the nucleon at an arbitrary scale requires 
a scheme which consistently interpolates between the two limits of VFNS and FFNS. A 
number of authors have pursued the construction of such schemes, both for unpolarized 
il] . ^8|pO| , |2T|| and polarized |22| scattering. 



With the exception of Ref. |2(|, where the splitting functions were modified for quark 



mass effects, the essence of the interpolating schemes is to incorporate a mass dependence, as 
a function of m 2 /Q 2 , in the coefficient functions which would enable Eq.([l|) to be recovered 
when Q 2 3> m 2 , and Eq.(^) when Q 2 ~ m 2 . By keeping the evolution equations unmodified, 
one ensures that the parton distributions (for the 3 light quarks, for the gluon, and for the 
newly introduced charm distribution) are defined consistently in the MS scheme. 

Until recently, quark mass corrections were only included in the gluon coefficient function, 
H g , while the quark coefficient function was taken to have the form derived for massless 
quarks [jT^jnfl . In fact, it is quite natural to think that the gluon sector, through the P cg 



splitting function, would be responsible for the generation and evolution of a heavy quark 
distribution. Furthermore, any quark mass correction (other than that appearing in the 
new, massive, scaling variable) would appear in the 0(a s ) coefficient functions only. In 



Ref. [15], Kretzer and Schienbein also calculated the quark mass corrections to the quark 
coefficient function. They found that these were small in absolute terms, but of a similar 
order of magnitude as the total in some regions of x, and for some choices of factorization 
scale. There are at least two reasons, therefore, why one needs to include these effects in 
the present analysis. Firstly, with the introduction of intrinsic charm the quark sector may 
become relatively more important than the gluon sector in the region of x where intrinsic 



charm is non-zero; and secondly, the analysis of Ref. |15| focused on quark mass corrections 
to the quark coefficient functions at relatively small values of x, while here we also need to 
examine the size of these corrections at large x. 

Of course, the full massive coefficient functions introduce extra mass logarithms in the 
quark and gluon sectors, resulting in double counting when using the evolution equations 
for the charm and gluon densities. This is reminiscent of a scenario where logarithms have 
not yet been resummed by the renormalization group equations. The problem can be cir- 
cumvented easily enough by requiring that in the limit Q 2 — > oo the massless coefficient 
functions are recovered after introducing subtraction terms H^ u g h : 

lim (H q (z,Q 2 ,m 2 c )-H™ b (z,Q 2 ,m 2 c )) = C q (z,Q 2 ) } (7a) 
lim (H g (z,Q 2 ,m 2 c )-H™ h (z,Q 2 ,m 2 c ))=C g (z,Q 2 ), (7b) 

where the factorization scale is chosen to be /i 2 = Q 2 for both the interpolating scheme and 
the VFNS (a choice to which we restrict ourselves throughout this analysis). The function 
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H q (z,Q 2 ,ml) in Eq.(|7a|) is the full massive quark coefficient function at 0(a s ), calculated 
in Ref. fl5[]. The subtraction terms are given by: 



* Z7T 



I±f!(i n 3!_l_21n(l-z) 
1 — z \ mi 



(8a) 



a s (Q 2 )^ , Q_ 2 

,2 



H^ b (z,Q\mi) = ^^T F ln^ [l-2z + 2z i ). (8b) 



Finally, one has: 

F 2 c (x, Q 2 ) = e 2 /' (fr 9 («, Q 2 , m 2 c ) - H™ b (z, Q\ m 2 )) * (c(£/*, Q 2 ) + eft/*, Q 2 )) 
+ e 2 / dz H g (z,Q 2 ,m 2 c ) -g(x/z,Q 2 ) 

J x Z 

- e 2 /' efe tf g sub (^,Q 2 ,m 2 ) -g{£/z,Q 2 ), (9) 
where the scaling variable ^ includes quark mass corrections, 

(10) 

The charm and gluon distributions in Eq.(|9]) are determined from evolution in exactly the 
same way as in the massless case described in Section II. 

Figure |2] shows F| evaluated according to the three schemes described above at x = 0.05 
and 0.2, for a Q 2 range relevant to the EMC data. In addition, for the interpolating scheme 
we also show the effect of neglecting the quark mass dependence in the quark coefficient 
function, so that H q -H^ nh is replaced by C q in Eq.(^|). For the parton distributions the 
GRV parameterizations ]13| are used. Note, however, that in the GRV fit a charm quark 




distribution is never introduced explicitly, rather F 2 C is always calculated via the FFNS. 
The charm density is generated from the GRV distributions by evolving with 3 flavors from 
/i 2 = 0.4 GeV 2 to m 2 in next-to-leading order, then from m 2 to Q 2 with 4 flavors according 
to the VFNS. For the FFNS calculation, the gluon distribution is evolved in leading order 
from fi 2 = 0.26 GeV 2 to m 2 c with 3 flavors. (Note that the choice of /i 2 in Figure |2| is only 
for the purpose of comparison with the other methods, which introduce charm at m 2 . In 
the final calculations the scale fi 2 = 4m 2 will be used for the FFNS, which is also the value 
used in Ref. 0.) 

At small x the effect of the mass-corrected quark coefficient functions on F 2 C turns out to 
be negligible, and only slight at larger x. As Q 2 becomes large, one can see in Figure |2] how 
the VFNS and interpolating schemes converge. Even at small Q 2 the difference between 
these is not large. On the other hand, while the FFNS provides a good approximation to 
the interpolating scheme for Q 2 < 30 GeV 2 , it dramatically overestimates the full result at 
larger Q 2 , especially at large x. Since this is the region where most of the relevant EMC 
data which we analyse lie, clearly we need to use the full interpolating scheme in order to 
draw reliable conclusions from our analysis. 

Another source of uncertainty in the calculation of F 2 C comes from the gluon distribution 
at large x, which at present is not very well constrained. To cover the full range of allowed 
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gluon distributions, we use the maximum and minimum gluon distributions from the MRST 
12 1 parameterizations in addition to that of GRV [|U|. The resulting F 2 C for the different 



next-to-leading order glue is shown in Figure [3], where evolution was again performed as 
described in the massless evolution section The data at small x in Figure |3] are taken from 
the ZEUS Collaboration |J , while the large-x data are from the earlier EMC experiment [|3] . 
At small x all of the parameterizations fit the data very well. At large x the maximum-gluon 
MRST and GRV fits also provide good descriptions of the data, with the exception of the 
two points at Q 2 = 60 GeV 2 . The last data point at x = 0.44 is not shown in the fits of Refs. 
1^.0], even though this point appears to confirm the trend indicated by the x = 0.24 point 



to lie somewhat above the perturbative QCD calculation. Neglecting the large-x, Q 2 = 60 
GeV 2 points, one would conclude that perturbative QCD fits the F 2 C data very well, without 
any need for additional non-perturbative contributions |12|,|13j. On the other hand, taking 



these points seriously has led several authors |||6| to conclude that the large- x EMC data 
provide evidence for an intrinsic charm component of the nucleon. 

In the next Section we shall study the large- a; EMC data more carefully, with the aim 
of ascertaining whether these can be understood perturbatively, or whether they can indeed 
be interpreted as suggesting that a perturbative QCD treatment alone is incomplete. 

V. NON-PERTURBATIVE CHARM 

The apparent discrepancy between some of the large-x F% data and predictions based 
solely on perturbative QCD has prompted several authors [j§-|6| to take seriously the possi- 
bility that an additional, non-perturbative, component of _F 2 C may be necessary to account 
for the data over the full range of x and Q 2 . In this Section we discuss how non-perturbative 
charm may affect _F 2 , particularly at large x, and how the intrinsic contributions can be 
included on the same footing as the perturbative effects. 

In earlier analyses intrinsic charm distributions have simply been added to the 
perturbatively generated F 2 , 

F<(x,Q 2 ) = F c 2 (pert) (x,Q 2 ) + F 2 C (IC) (*,Q 2 ), (11) 

where the perturbative contribution, F 2 ( - pert - ) ) is given by Eq.(|4]) (with higher order correc- 
tions), while the intrinsic charm contribution, F 2 ^ IC \ in its simplest form is: 

F 2 c(IC) = e 2 x(c IC + c IC ). (12) 



In practice, 0(a s ) contributions to Eq.(12) are fully implemented in this analysis (the rele- 
vant expressions are given in Ref. |J). 

Within the interpolating scheme of Section IV, the most natural way to implement 
intrinsic charm in F 2 C is to modify the boundary condition for the charm quark distribution. 
Instead of c(x, /i 2 ) = c(x, /i 2 ) = 0, one now has non-zero distributions at the scale /i 2 = m 2 . 
The physical reasoning for this is that if there are non-perturbative processes producing 
charm in the nucleon, this charm can be resolved (brought on its mass shell) only when the 
system has sufficient energy. In the MS scheme, the scale at which the number of flavors 
changes from 3 to 4 is fj, 2 = m 2 , so that regardless of the dynamical origin of the charm, there 
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will be enough energy in the system to open a new active flavor channel for Q 2 > m 2 . With 
this in mind, we next discuss several non-perturbative models which attempt to describe the 
generation of intrinsic charm in the nucleon. 



A. Five-Quark Component of the Nucleon (IC1) 



Based on the initial observation |23| that the charm production cross section in hadronic 
collisions was larger than that predicted in leading order perturbative QCD, Brodsky et 
al. M suggested that the discrepancy could be resolved by introducing an intrinsic, non- 
perturbative, charm component in the nucleon wave function. In this model, which we 
shall refer to as "IC1", the nucleon is assumed to contain, in addition to the lowest energy 
three-quark Fock state, a more complicated, five-quark configuration on the light-cone: 

\p) = co\uud) + ci\uudcc), (13) 

where Cq (cf) is the three (five)-quark probability. In order to explain the original data 



23] , the normalization of the latter was chosen to be 1% 0]. Assuming the five-quark wave 
function to be inversely proportional to the light-cone energy difference between the nucleon 
ground state and the five-quark excited state, one finds that the x-dependence of the c quark 
distribution is given by [|J: 

c IG1 (x) = 6x 2 ((1 -x)(l + 10x + x 2 ) - 6x(l +x)logl/x) . (14) 



The anticharm distribution has the same shape as the charm distribution in this model, 



c lcl (x) = c lcl (x) 



Because the intrinsic charm in this model is assumed to be generated through gg — > 
cc processes, with each gluon originating from different valence quarks, the cc probability 
scales like a 2 (m 2 )/m 2 relative to the perturbative component PJ24"! . This contribution was 
therefore interpreted in Refs. PJ24| as a higher- twist effect. Our philosophy here is that 
since this is generated non-perturbatively, the resulting non-perturbative intrinsic charm 
distribution calculated at m 2 should be evolved as leading twist, on the same footing as 
the perturbative contribution. Since there is only one kind of charm quark, irrespective 
of its origin, QCD corrections should affect both the perturbative and non-perturbative 
distributions identically. 



B. Meson Cloud Model (IC2) 



An alternative to the five-quark intrinsic charm model was considered in Refs. [25-27 



in which the charm sea was assumed to arise from the quantum fluctuation of the nucleon 
to a virtual D° + configuration. In the following we shall refer to this model as "IC2" . 



The nucleon charm radius |25| and the charm quark distribution ]26| were both estimated 
in this framework. Furthermore, the effects of hard charm distributions on large- x HERA 
cross sections, and in particular on the so-called HERA anomaly H28|j29|| , were studied in 
Ref. IE7I- 



The meson cloud model for the long-range structure of the nucleon has been used exten- 
sively to describe various flavor symmetry breaking phenomena observed in deep-inelastic 
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scattering and related experiments. It offers a natural explanation of the d excess in the 
proton over u [§l],[n]] in terms of a pion cloud, which itself is a necessary ingredient of the 
nucleon by chiral symmetry. It also provides an intuitive framework to study the strangeness 
content of the nucleon, through the presence of the kaon cloud |32|| . Whether the same phi- 
losophy can be justified for a cloud of heavy charm mesons and baryons around the nucleon 
is rather more questionable given the large mass of the fluctuation. Nevertheless, to a crude 
approximation, one may take the meson cloud framework as indicative of the possible order 
of magnitude and shape of the non-perturbative charm distribution. Furthermore, a natural 
prediction of this model is that the c and c distributions are not symmetric. 

In the meson cloud model, the distribution of charm and anticharm quarks in the nucleon 
at some low hadronic scale can be approximated by |27f : 



IC2. 



c 



c IC2 (x) « f D/N (x), (15b) 



where 



f - , (x) - — r dk 2 g2{kl > x) ( kj + [M Ac -(l-x)M]2 \ 

is the light-cone distribution of D° mesons in the nucleon, and /a c /7v( x ) = fjjm(l — x) 
is the corresponding distribution of A+ baryons. In Eq.(|l6|) the function g describes the 
extended nature of the DA C N vertex, with the momentum dependence parameterized by 
g 2 (kj_, x) = g 2 (A 2 + M 2 )/(A 2 + s), where s is the DA C center of mass energy squared and 
go the DA C N coupling constant at the pole, s = M 2 . As a first approximation, one might 



take go to be of the same order of magnitude as the ttNN coupling constant. In Ref. |33 
this coupling constant was estimated within a QCD sum rule calculation. 



C. Intrinsic Charm Distributions 

The c and c distributions in the intrinsic charm models IC1 and IC2 are shown in Figure |], 
each normalized to a common value of 1%. For the IC2 model this corresponds to a cutoff 
A w 2.2 GeV (for a probability of 0.5%, for example, one would need A w 1.7 GeV). Quite 
interestingly, the shapes of the c quark distributions are quite similar in the two models, 
with xc peaking at around x ~ 0.3. However, because the IC2 model gives a significantly 
harder c distribution, while IC1 implies that c and c are equal, the resulting structure 
function, F£, will be somewhat harder in the IC2 model. For comparison, a typical (soft) 
perturbatively-generated charm distribution is also shown in Figure |], evaluated from the 
MRST parameterization |12| (with the maximal gluon) at Q 2 = 5 GeV 2 . 

The effects of the modified boundary conditions incorporating non-zero intrinsic charm 
distributions are shown in Figure for the GRV parameterization [H| and for 1% intrinsic 



charm normalizations, at Q 2 = 25, 45 and 60 GeV 2 . The data at 60 GeV 2 are well fitted with 
a 1% IC1 component, although with the IC2 model one slightly overestimates the x = 0.44 
data point, due to its very hard c distribution. At lower Q 2 values, however, the addition 
of a 1% intrinsic charm component, from either model, overestimates the large-x points. 
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This finding is essentially independent of the parton distribution functions employed, as 
Figure |6| illustrates for the MRST distributions [Q. From this one can conclude that with a 
1% intrinsic charm component one cannot simultaneously resolve the large-x discrepancy for 
the large-Q 2 data, and maintain a satisfactory fit to the data at lower Q 2 . One could imagine 
introducing a Q 2 -dependent normalization for the intrinsic charm to improve the fits to the 
25 and 45 GeV 2 data, although this solution would not be well grounded theoretically. 

To compare with the procedure for incorporating intrinsic charm adopted in the earlier 
analysis in Ref. ||, we show in Figure ^ the F% obtained from the FFNS through Eq. (|TT|) , 
and the 0(a s ) corrections to the intrinsic charm component, F% ^ — see Eq.(57) of Ref. 
0. (Note that although the FFNS curves appear to lie slightly below the data for Q 2 = 25 
and 45 GeV 2 , inclusion of 0(a 2 ) corrections leads to a slight improvement for the lower Q 2 
data, without much effect on the data at 60 GeV 2 []13[ .) The results are qualitatively similar 
to those obtained from the full interpolating scheme, namely the data at different Q 2 seem 
to require different functional forms and normalizations for the intrinsic charm. 

By varying the amount of intrinsic charm that the data can accommodate, one can 
attempt to improve the overall fit at all values of Q 2 . Using the minimum gluon fit of the 
MRST parameterization one can determine the maximum intrinsic charm allowed by the 
data for both the IC1 and IC2 models. Slightly better fits can be obtained with a 0.75% IC1 
charm distribution, and with a 0.4% IC2 distribution, as indicated in Figure ^ for Q 2 = 45 
and 60 GeV 2 . (The quality of the fit at Q 2 = 25 GeV 2 is similar to that for the 45 GeV 2 
data.) However, even given the different shapes of the charm and anticharm distributions 
in the IC1 and IC2 models, it is still quite difficult to obtain a quality fit to data at all 
Q 2 values. For the Hoffman and Moore procedure || with the FFNS, one finds similar 
difficulties in reconciling the Q 2 = 60 GeV 2 data with those at lower Q 2 , even as a function 
of the intrinsic charm model, and normalization (c.f. Figure ||), although there may be a 
slight preference for IC2 with 0.4% normalization. 

In view of the difficulties in obtaining a simultaneous fit to the F£ data at all measured 
values of x and Q 2 with either the perturbative-only, or the intrinsic charm scenarios, it 
seems imperative that the data, particularly those at large x and Q 2 , be confirmed by 
independent measurements. 



VI. CONCLUSION 

Unraveling the rich structure of the nucleon sea is an ongoing endeavor which has taught 
us a number of important and sometimes surprising lessons in recent years. Far from being 
a homogeneous background in which the valence quarks reside, the sea has proved to exhibit 
considerable intrinsic structure of its own. Perhaps the most conspicuous structure is that 
in the light quark sea of the proton, where a number of experiments have now confirmed 
beyond any doubt a significant excess of d over u quarks ||30|| . This in turn has created an 



environment in which the importance of non-perturbative effects in the nucleon is appreciated 
to a far greater extent, even when discussing structure at deep-inelastic energy scales ]3j] . 

More speculative, and less constrained experimentally, are suggestions that the proton sea 
for heavier flavors might also exhibit characteristics which cannot be attributed to perturba- 
tive QCD mechanisms alone p7|j32]| . A prime example would be the presence of asymmetric 



sea quark and antiquark distributions, which have been searched for in the strange sector in 
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both deep-inelastic neutrino and antineutrino scattering, as well as in electromagnetic form 
factors at low energies |p4 |. Indeed, there is no symmetry principle in QCD which would 



prevent sea quarks and antiquarks having different momentum distributions, just as there 
is no symmetry requiring the d and u sea to be equivalent. 

More challenging still is to identify possible non-perturbative effects in the charm sector, 
especially since the magnitude of the total charm component of non-charm hadrons is tiny. 
Nevertheless, a number of pioneering studies of various charm production reactions 
have left us with lingering doubts as to whether perturbative QCD is the full story behind 
these processes. On the other hand, the evidence doesn't appear to be conclusive enough 
to warrant introduction of intrinsic charm distributions in global data parameterizations 



12| . ^3|j35| , which to date have generated the charm structure function purely perturbatively. 

The present analysis has been an attempt to resolve the issue of whether the charm elec- 
troproduction data do indeed support the existence of non-perturbative components of the 
nucleon wave function, or whether they can be understood within conventional perturbative 
dynamics. To this effect we have used the latest available technology to describe charm 
production over the entire range of x and Q 2 accessible to experiment. Our approach con- 
sistently interpolates between the two asymptotic regions of massless evolution at large Q 2 
and the photon-gluon fusion process at low Q 2 , and includes quark and target mass effects 
and corrections for mass thresholds. 

To a certain extent our findings confirm the existing state of affairs, in which some of the 
data show no evidence at all for intrinsic charm, while other data cannot be fitted without 
additional non-perturbative input. Even within the rather different models of intrinsic charm 
considered here, with varying normalizations, it seems difficult to simultaneously fit the 
entire data set in terms of a single intrinsic charm scenario — although there may be a 
slight preference for intrinsic charm in model IC2 at a level of about 0.4%. The clearest 
conclusion that one can draw from this is that more quality data are urgently needed to settle 
the issue. In particular, while the small- X domain seems to be relatively well under control, 
the large- X region, where cross sections are small and measurements more difficult, must be 
the focus of future experimental effort if the non-perturbative structure of the nucleon sea 
is to be resolved. 
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FIGURES 
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FIG. 2. Charm structure function as a function of Q 2 for x = 0.05 and 0.2, evaluated according 
to the various schemes discussed in the text. "IS" refers to the interpolating scheme of Eq.(||), while 
"IS (massive g)" denotes the interpolating scheme but with massless quark coefficient functions. 
The GRV parameterization [13| for the gluon is used. 
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FIG. 3. Charm structure function calculated within the interpolating scheme for different gluon 
distributions, GRV |jl3| (dotted) and MRST |[^] with minimum (dashed) and maximum (solid) 
gluons. The data at small x are from the ZEUS Collaboration Q, while the large- x data are from 
the EMC ||. For clarity the small- x curves have been scaled by a factor 10 (0.1) for Q 2 = 25 (60) 
GeV 2 , and the large- x curves by a factor 100 (0.01). 
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FIG. 4. Charm quark distributions from the intrinsic charm models IC1 |27j (solid) and IC2 
(dashed), both normalized to 1%, and from the MRST parameterization |l2j (with maximal 
duon) at Q 2 = 5 GeV 2 (dotted). 
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FIG. 5. The charm structure function calculated through the interpolating scheme ("IS") in 
Eq.(H|) for intrinsic charm distributions from models IC1 and IC2, normalized to 1%. The GRV 
parameterizations for the gluon and light quark densities are used. 
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FIG. 8. As in Figure ^, but with different normalizations for the IC1 and IC2 model distribu- 
tions, and with the MRST parameterization [O] with the minimum gluon. 
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